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Abstract. We are concerned with the critical threshold phenomena in the Restricted Euler 
(RE) equations. Using the spectral and trace dynamics we identify the critical thresholds for 
3D and the 4D restricted Euler equations. It is well known that the 3D RE solutions blow 
up. Projected on the 3-sphere, the set of initial eigenvalues which give rise to bounded stable 
solutions is reduced to a single point, which confirms that 3D RE blowup is generic. In contrast, 
we identify a surprisingly rich set of the initial spectrum on the 4-sphere which yields global 
smooth solutions; thus, 4D regularity is generic. 

To Katepalli Sreenivasan on his 60th birthday, with friendship and admiration 

1. Restricted Euler Equations and Spectral Dynamics 

We are concerned with the questions of global regularity vs. finite-time breakdown of Eulerian 
flows governed by 

d t u + u-V x u = F, t>0. 

Here u is the velocity field, u := (u±, 112, ■ ■ ■ , u n ) T : R 1+n 1— ► R n , and its global behavior is 
dictated by the different models of the forcing F = F(u, Vu, . . .). For forcing involving viscosity 
and pressure, we meet the well known Navier-Stokes (NS) equations, 

(1.1) d t u + u ■ V x u = vAu - Vp, x E M n , t>0, 

augmented with the incompressibility condition, V • u = and subject to prescribed initial 
conditions u(x, 0) = uq(x). In many applications, v > is sufficiently small so that one can 
anticipate the behavior of slightly viscous NS solutions to be described by the Euler equations 
with v = in (jl.ip . at least for flows occupying the whole space so that the important effects of 
boundary layers can be ignored. 

The velocity gradient of the incompressible Euler equations, M := V x u then solves 

(1.2) dtM + u- V X M + M 2 = -(V® V)p. 

Taking the trace of (jl.2p while noting that M is trace-free, trM = V • u = 0, one finds trM 2 = 
— Ap which dictates the pressure as p = — A~ l {trM 2 ). The second term in (11. 2p therefore 
amounts to the n x n time-dependent matrix 

(V (8) V)A' 1 {trM 2 ) = R[trM 2 }. 

Here R[w] denotes the so called Riesz matrix — an n x n matrix whose entries, (R[w])jk ■= 
RjRk(w), involve the Riesz transforms Rj, Rj = — (— A)~ l l 2 dj, i.e., 

R[w]:={RjR k (w)}l k=1 , R^W(0 = for 1 < j < n. 
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This furnishes an equivalent, self-contained formulation of Euler equations, expressed in terms 
of the velocity gradient M = V x u, which is governed by, 

(1.3) d t M + u- V X M + M 2 = R[trM 2 }, M = V x u, 

and subject to the trace- free initial data, M(-,0) = Mq. Observe that the invariance of in- 
compressibility is already taken into account in (|1.3|) since trM 2 = trR[trM 2 ] implies (dt + u ■ 
V x )trM = and hence trM = trM = 0. 

It is the global nature of the Riesz matrix, R[trM 2 ], which makes the issue of regularity for 
Euler and NS equations such an intricate question to solve, both analytically and numerically, 
[6]. Various simplifications to this pressure Hessian, R[trM 2 ] = — V (g) Vp were sought, e.g., 
[TP] [TBI 01 El [2] . In this paper we focus our attention on the so called restricted Euler equations, 
proposed in (9[ [16] as a localized alternative of the full Euler equations (jl.3p . By the definition 
of the Riesz matrix, one has 

R[trM 2 } = V ® VA _1 [trM 2 ] = V tg> V / K(x - y){trM 2 )(y)dy, 

where the kernel K(-) is given by 

1 



K(x) 



■ln\x\ n=2, 

Z7T 



1 



n > 2, 



{(2 - n)u n \x\ n ~ 2 

with oj n denoting the surface area of the unit sphere in n-dimensions. A direct computation 
yields 

n a\ a a is + trM 2 f \x - y\ 2 5 jk - n(xj - yj)(x k - y k ) A 

(1.4) djd k K*trM' = 5 jk + / J - . J — — % trM\y)dy. 

n J R n u n \x-y\ n+2 

Ignoring the singular integrals on the right of (jl.4p , we are left with the local part of the Riesz 
matrix R[trM 2 ], given by trM 2 I nxn /n. We use this local term to approximate the pressure 
Hessian in (jl.3p . The resulting Restricted Euler (RE) equations amount to 

(1.5) d t M + u ■ V X M + M 2 = —trM 2 I nxn . 

n 

This is a matrix Ricatti equation for the n x n matrix M, which should mimic the dynamics 
of the velocity gradient, Vu in the full Euler equations. We observe that as in the full Euler 
equations, incompressibility is maintained in the restricted model, since trM 2 = tr[trM 2 I nxn /n] 
implies (d t + u- V x )trM = and hence trM = trM = 0. The 3D RE (fT3[) has attracted great 
attention since it was first introduced in j9[ [16] as a local approximation to the full 3D Euler 
equations. It can be used to understand the local topology of the Euler dynamics and to capture 
certain statistical features of physical turbulent flows, consult [Tl [51 [T6]. 

What about the global regularity of the RE equations (|1.5|) ? the finite time breakdown of the 
three-dimensional restricted model goes back to the original work of Viellefosse |16| . In [TTJ we 
have shown that the 3D RE solutions break down at a finite time for all initial configurations 
Mo, except for the special case when Mq has three real eigenvalues, 

Ai(0) < A 2 (0) < A 3 (0), {Aj-(O) = \i(M )} 3 . =1 , 

which are aligned along the ray (— r, — r, 2r), r £ M + . Thus, the finite time break down of the 
3D RE equations is generic. 

In this paper we shall identify and compare between the restricted Euler equations in 3D and 
4D case, respectively. To this end, we consider a bounded, divergence-free, smooth vector field 
u : M. n x [0, T] — > W 1 . Let x = x(a, t) denote an orbit associated to the flow by 

dx 

— = u(x, t), 0<t<T, x(a, 0) = a G R n . 
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Then along this orbit, the velocity gradient tensor of the restricted Euler equations (|1.5j) satisfies 

—M + M 2 = I nxn , —:=dt + u-V x . 

at n at 

By the spectral dynamics lemma 3.1 in |llj . the corresponding eigenvalues of M satisfy 

(1.6) ^ + ^ = -£4 < = l,-,n- 

i=i 

This is a closed system for A(t) = (Ai(i), ^(t), . . . , X n (t)), which serves as a simple approxima- 
tion for the evolution of the velocity gradient field. 

For arbitrary n > 3, we use the spectral dynamics of M in order to show the existence of 
a large set of initial configurations leading to finite time breakdown of (|1.6|) . generalizing the 
previous result of [16]. The finite time breakdown of the n-dimensional RE equations (and the 
precise topology of the breakdown) was established in [TT] after we identified a set of [n/2] + 1 
global spectral invariants, interesting for their own sake. Yet, this does not exclude the possible 
existence of other generic sets of initial data, for which global smooth solutions exist. The 
distinction between these two sets of initial conditions is identified by the so-called critical 
threshold surfaces in configuration space: finite time breakdown occurs for super-critical initial 
data on "one side" of the such critical threshold, while the set of sub-critical initial data on the 
"other side" of the threshold yields global smooth solutions. 

An interesting question therefore arises, namely, whether there exists a critical threshold for 
the 4D restricted Euler equation. This remarkable critical threshold phenomena was identified 
in [17] for a class of essentially ID Euler-Poisson equations, and in [12] for a ID convolution 
model for nonlinear conservation laws. The 2D critical threshold phenomena has been recently 
confirmed for a restricted Euler-Poisson system |13j and a rotating Euler equation [141 E] • I n an 
these cases, we identified large, generic sets of sub-critical initial data, which evolve to global 
smooth solutions. This is in contrast to the generic scenario of finite-time blows up in the 3D 
RE equations^. 

In this paper we identify the exact critical thresholds for the 4D restricted Eulerian (RE) 
equations and we conclude with the surprising result that in the 4D case, the RE equations 
admit a large, generic set of sub-critical initial data which give rise to global smooth solutions. 

A summary of our results is outlined below. We say that Ao G M n is sub-critical if there exists 
a global solution in time of (jl.6p . subject to initial conditions, A(0) = Ao. A first observation 
rests on the obvious symmetries of f)1.6|) . 



Lemma 1.1. If A is sub-critical then so is rA, Vr > 0. Moreover, A a = {A^q), Vcr G vr n } is 
also subcritical. 

For the proof we note that if A(t) is the global solution corresponding to Ao, then rA(ri), r > 
is the global solution corresponding to rAo- Also, equations (II. 6|) remain invariant under 
arbitrary permutation a which amounts to reordering, exchanging the Aj-equation with X a (jy 
equation. It follows that the set of sub-critical initial data consists of rays, and therefore, it is 
enough to consider the projection of this set on the unit sphere. In fact, we can restrict attention 
to an orthant of any convex set containing the origin. In this context we have 

Theorem 1.1. Solutions to hi. 0^1 with n = 3 remain bounded for all time if and only if the 
initial data Ao := (Aio, A20) A30) lie in the following set 

r{(-l,-l,2) CT }. 



^We should emphasize that generic sub-critical data are not limited to a perturbative statement of global 
existence for initial data in the local neighborhood of certain "preferred configurations". Instead, the precise 
notion of "generic" sub-critical sets, quantified below and the references mentioned above, makes it clear the 
critical threshold phenomena we seek describes a global scenario in configuration space. 
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Restricted to one orthant of the unit sphere, we thus find that the 3D RE equations admit only 
one sub-critical point. In this sense, the finite-time breakdown of 3D RE is generic. This result 
was already obtained in pTJ by spectral dynamics analysis. In §4 we present an alternative, 
equivalent argument based on trace dynamics of tr(M k ), 1 < k < n, which paves the way for 
identifying our 4D critical threshold surface in §5. 

In contrast to this generic 3D finite-time breakdown, the 4D RE equations admit a large class 
of global smooth solutions. Our 4D results are summarized below. 

Theorem 1.2. Solutions to $1.6}) with n = 4 remain bounded for all time if and only if the 
initial data Aq := (Aio, A20J A30, A40) with Y2j=i Ajo = 0, up to a permutation, lie in one of the 
following sets 

{i}Two pairs of arbitrary complex eigenvalues, 

Aq e {(a + ib, a — ib, —a + ic, —a — ic), be 7^ 0}. 

{ii} One pair of complex eigenvalue with two equal real eigenvalues 

Aq G {(a + ib,a — ib, —a, —a), b / 0}. 

{Hi} Real eigenvalues, 

Aq £ {(a + b, a — b, —a, —a), b G [—2a, 2a], a > 0}. 

Expressed in terms of traces := tr(M k ), these initial configurations form a "large" sub- 
critical set which can be realized by its projection on the surface E, 

4 

E := {A I 4m 4 - 2ml ~ 2m 2 + 3 = 0, m 1 = 0}, m k : = ^ A*. 

k=l 

Theorem 1.3. Solutions to U.6\) with n = 4 remain bounded for all time if and only if there 
exists a r > such that the initial data Aq := (Aio, A20) A30, A40) lie in the following set 

rA 6 Sfl [{|m 3 | < 1.5(1 - m 2 ),m 2 < 1} U {m 3 = 1.5(?n 2 - 1), m 2 > 1}] . 

The set stated in Theorem 11.31 is non-trivial; in fact, it contains non-zero neighborhoods. We 
note that a recent study in [15] on 4D incompressible Navier-Stokes equations suggests that the 
energy transfer in 4D is indeed more efficient than in 3D. 

After this introduction of restricted Euler equations and the associated spectral dynamics, we 
identify the 4D sub-critical initial configurations in terms of eigenvalues in §2. An alternative 
formulation of the spectral dynamics — called trace dynamics is derived in §3. Based on the 
trace dynamics we identify the critical thresholds for 3D case in §4 and the 4D model in §5. 
Finally in the appendix we establish the correspondence between the sub-critical eigenvalues 
and the sub-critical set for initial traces. 

2. 4D SPECTRAL DYNAMICS 

Let A = Aj solve the restricted Euler equation 

(2.1) ^ + A? = ^X>f, i = l-.4. 

Two independent global invariants obtained in [JTJ are 

(2.2a) (Ai - A 2 )(A 3 - A 4 ) = (Aio - A 20 )(A 3 o - A 40 ) 

and 



(2.2b) 



(Ai — A3XA2 — A 4 ) — (Aio — A30XA20 — A40) 
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We now prove Theorem 11.21 based on these global invariants. In view of the incompressibility 
invariant Y2j=i -\? = 0' we can express the remaining three spectral degrees of freedom as 
A = (a + b, a — b, —a + c, —a — c) T , where a is real, a£R, and b or c are either real b, c £ K or 
purely imaginary, b, c £ iK. The two global invariants (|2.2p now read 

(2.3a) be = &0C0 

and 

(2.3b) 4o 2 - b 2 - c 2 = 4a 2 , - b 2 Q - c 2 ,. 

The spectral dynamics (|2.1|) amounts to the 3x3 closed system, 

(2.4a) la = --^ + -c 2 , 

y 1 dt 2 2' 

(2.4b) j t b = -lab, 

(2.4c) — c = -2ac, 

at 

subject to initial data (ao, &0) c o)- Observe that both 6 = and c = are global invariants, thus 
the only equilibrium points of (|2.4p when &o c o 7^ lie along the curves (0, b*,c*), b* = ±c*. 
From (|2.4bp and (|2.4c|) it is clear that if either bo or cq are purely imaginary then they remain 
so for all time. Thus, we need to discuss three cases in order. 

{i} Two pairs of complex eigenvalues, ao ± ibo and — ao ± icq with boco ^ 0. Setting (a, b, c) i— ► 
(a, i&, zc) in (|2.3bj) we obtain the global invariant 

4a 2 + b 2 + c 2 = 4a 2 , + 6q + Co- 
in this case, all trajectories remain bounded for all time. 

{ii} One pair of complex eigenvalues, aoiibo, bo / and — ao±co- Setting (a, b, c) i— > (a, z6, c) 
in (|2.3bp . then the global invariant (|2.3b|) becomes 

4a 2 + 6 2 -c 2 =Aa 2 + bl-cl. 

We distinguish between two cases. If cq = 0, then by (|2.4c[) c(t) = 0, and the reduced global 
invariant, 4a 2 + b 2 = 4a 2 , + &q, implies that both a and b remain bounded for all t > 0. If co 7^ 0, 
then the equation (|2.4ap becomes 



dt 2 V & 2 

this shows that no finite equilibrium point of the system (|2.4a|) - (j2.4c|) is stable, which excludes 
the possibility of a globally bounded solution when 60C0 7^ 0. 

{hi} Two real eigenvalues, ao ± 60 an d — ao i c o- Again, we distinguish between two cases. 
Assume that two initial eigenvalues coincide, say cq = (if 60 = 0, we end up with a similar 
scenario which amounts to a permutation of the cq = case). Then c(t) = and the remaining 
(a, b) satisfy the reduced 2x2 system 

d ( a\ _( -b 2 /2 



dt\b J \ -2ab 

with the corresponding global invariant 4a 2 — b 2 = \a\ — b\. Now, since = —b 2 /2 < 0, it 
follows that a > is decreasing while 6 must approach the stable equilibrium points (a* > 0, 0) 
along the positive a-axis, as ~^b = —2ab has the opposite sign of b. Thus, trajectories remain 
bounded in the invariant sector |&o| < 2ao- 

Finally, assume no pair of initial eigenvalues coincide, 6o c o 7^ 0. Then, since the global 
invariants (|2.3ap and (|2.3bp are not compact, the only possible bounded solutions are those 



6 HAILIANG LIU, EITAN TADMOR, AND DONGMING WEI 

converging to the equilibrium points (0, ±c*,c*). But when substituted into both (|2.3ap and 
(l2l3b"l) . this impl ICS 

4a 2 , = (b ± c ) 2 , 

which is satisfied only if at least one pair of initial eigenvalues coincide, i.e. ao ± bo = — a o ± cq. 
We conclude that for real eigenvalues, only those of the form (ao + &o>«o — bo, —ao, — ao) with 
I fro I < 2ao,«o > lead to global bounded solution. 

3. Trace dynamics 

This section is devoted to an alternative formulation of the spectral dynamics in terms of real 
quantities m k := Y^j=i > k = 1, - ■ ■ ,n, where A = Aj solves the restricted Euler equation 

(3.1) ^ + ^£A 2 , < = l-«- 

j'=i 

This is motivated by the trace dynamics originally studied in [16] for n = 3. The use of trace 
dynamics enables us to obtain an explicit description of the critical threshold surface for initial 
configurations. 

Here we seek an extension for the general n-dimensional setting, which is summarized in the 
following 

Lemma 3.1. [11] Consider the n- dimensional restricted Euler system \3.1\) subject to the in- 
compressibility condition mi := X^=i -\? = 0- Then the traces m k for k = 2, • • • ,n satisfy a 
closed dynamical system, see 113. ty) - [37b)) below, which governs the local topology of the restricted 
flow. 

Proof. Based on the spectral dynamics the evolution equation for each eigenvalue Aj can be 
written as 

d , ,2 1 

— Aj + Aj = -m 2 , i = l---n. 
dt n 

By multiplying kX^ -1 and summation over i we obtain 

d , k 

— m k + km k+ x = -m 2 m fc _i, k = 2---n. 
dt n 

Note that mi = we have 

(3.2) ^m 2 + 2m 3 = 0, 

dt 

(3.3) ^-m 3 + 3m 4 = -m 2 ., 

or n 

(3.4) 4: m n + ™m n+ i = m n _im 2 . 

To close the system, it remains to express m n+ i in terms of (mi, • • • , m n ). To this end we utilize 
the characteristic polynomial 

(3.5) X] + Ql X]~ l + ■■■ q^tXj +q n = 0, 
expressed in terms of the characteristic coefficients 

qi = -mi = 0, q 2 = -^m 2 , 03 = -m 3 /3, q 4 = -m 4 /4 + m|/8, 

Note that the q's can be expressed in terms of (mi, • • • , m n ). Using (|3.5p one may reduce m n+ i 
in (|3.4p to lower-order products. In fact, Y!ij=i(^j x P-^P j) gives 

(3.6) m n+ i + g 2 m n _i H h q n -im 2 = 0. 
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Substitution into (13, 4p yields the closed system we sought for. □ 

We demonstrate the above procedure by considering the two examples of 3D and 4D critical 
thresholds. 

4. 3D CRITICAL THRESHOLDS: FINITE TIME BLOWUP 

This section is devoted to the study of the 3D critical thresholds, see [161 HI- I n the three 
dimensional case one has 

1 " 1 

Qi = 0, q 2 = ~2 m 2, <?3 = Aj = --m 3) 

i=i 



hence 



A 3 rri2\i ma = 0, i = 1, 2, 3. 



1 . 1 



2 3 

Multiplying by Aj and taking the summation over i we find 

1 2 
TO 4 = 2 m 2- 

Thus a closed system is obtained, 

(4.1) 4 m 2 + 2m 3 = 0, 



dt 

(4.2) -|m 3 + = 

From p7T|) - (|3~2]) it follows that 



^[6m| - ml] = 6m 3 ^m 3 - 3m2^m 2 = 0, 



which yields a global invariant 



6m| — mj = Const. 



We consider the phase plane (777-2, 711,3), except for the separatrix 67713 = 777-3, an other solutions 
would not approach the origin. The phase plane is divided into two parts by this separatrix. The 
nonlinearity ensures that trajectories which do not pass the origin must lead to infinity at finite 
time. In fact for initial data from the region {(771-2,7773), 77i 2 > v^6 m^ 3 }, the corresponding 
trajectories will remain in this region since the system (|4.ip . (|4.2p is autonomous. Therefore 
T2D leads to 

d 1 3 /T7 4/3 



(4.3) — 7773 < — 2 V 36 in 

Since im^ = — |m|, vn^ (t) is always decreasing in time. Even for positive 7713(0), there exists 



a finite time T* such that 77i3(T*) < 0. The integration of (14. 3D over [T*,t) gives 

7113(7-) < 



3 1~ 3 
|^36(t-T*) +77i 3 (T*)- 1 / 3 

This shows that 7713 (t) — ► —00 when i approaches a time before 



T* + ^=(-77i 3 (T*))- 1 / 3 . 

3</36 



Finite time breakdown can be similarly justified for initial data lying in the region {(7112,1113), ui2 < 

4' 



■^6 m^ 3 }. These facts enable us to conclude the following 



8 HAILIANG LIU, EITAN TADMOR, AND DONGMING WEI 

Theorem 4.1. Consider the system fr4-l\ )- [4-W with initial data (7712(0), 7713(0)). The global 
bounded solution exists if and only if the initial data lie on the curve 



(7712,777,3) 



1 3/2 
7713 = 71 



We now turn to interpret this condition in terms of the eigenvalues. Set A = (Ai, A2, A3), the 
above critical stable set can be written as 

fc=i vu \fc=i / k=i 

The homogeneity of the above constraint in terms of eigenvalues implies that if A E 0, then 
rAefi Vr > 0. 

Without loss of generality we consider the restriction of Q, onto a ball Ylk=i ^1 = r2 ' denoted 
by f2(r). There are two cases to be considered: 

The initial eigenvalues contain complex components, say Ao = (a — bi, a + bi, c) for real 
a, b, c 6 M. The restricted set f2(\/6) is determined by 

c + 2a = 0, 2a 2 - 26 2 + c 2 = 6, 2a(a 2 - 3b 2 ) + c 3 = r 3 /^ = 6. 

Eliminating c we have 

6a 2 - 26 2 = 6, -6a(a 2 + b 2 ) = 6 4a 3 - 3a + 1 = (a + l)(2a - l) 2 = 0, 

which has real roots a S {—1,0.5,0.5}, from which no real b 7^ can be found. 

The only possible scenario is the real eigenvalue Ao = (a, 6, c) S M 3 . Restriction again on 
Vt we have 

a + b + c = 0, a 2 + 6 2 + c 2 = 6, a 3 + 6 3 + c 3 = r 3 /V6 = 6. 

Eliminating a, b we have c 3 — 3c— 2 = with real roots c E {2, —1,-1}. The symmetric property 
implies that a, 6 also lie in the set {2, —1, —1}. In short one has 

fi(V6) = {A|(-l, -1, 2), (-1, 2, -1), (2, -1, -1)}. 
This when combined with the above scaling property leads to the result stated in Theorem 1.1. 

5. 4D CRITICAL THRESHOLDS: GLOBAL REGULARITY 

In the four dimensional case one has 

1 1 777,4 

qi =0, q 2 = ~2 m 2> °3 = ~3 m 3, Qa = — ^ + 

Hence 

Af " ^A 2 - -m 3 A, t - + = 0, 7 = l---4. 
Multiplying by Aj and taking the summation we obtain 

1 1 5 

m 5 = -m 2 m 3 + -m 3 m 2 = -m 2 m 3 . 

2 3o 

Therefore the resulting closed system becomes 

d 

(5-1) —777,2 = -2m 3 , 

at 

d 3 9 

(5.2) —777,3 = ^"72 - 37774, 

^ ^ d 7 

(5.3) — m 4 = --77737712. 
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From (|5.ip and (15.31) it follows that 

d ( 7 2 \ 

— tua m 9 = 0, 

dt\ 12 2 J 

which gives a global invariant which gives the global invariant 

7 2 7 

u™ 2 = 77140 " n 1 



(5.4) m 4 -^i = ™ 40 -^!, 



Substitution of this into (|5.2p leads to 

^m 3 = - m | - 3(m 4 - Y^rnl). 

In order to ensure global bounded solution (excluding globally decreasing 7773) it is necessary to 
consider trajectories for which 

7 Z 2 

(5.5) m 4 (t) - — m|(i) = --, Z > 0. 
We thus have a closed system for (7772,7773) 

(5.6) = -2ms, ^rn 3 = -m\ + Z 2 

with a moving parameter Z determined by (|5.5p with i = 0. This system has two critical points 
(— Z, 0) and (Z, 0); it is easy to verify that as equilibrium points of system (15. 6p . (— Z, 0) is a spiral 
and (Z, 0) is a saddle for the corresponding linearized system. 

This structure suggests that part of separatrix' of this system may serve as the critical thresh- 
old. Note that 

— (3777,0 — 777,0 + 3Z 7772) = 6?773— 7773 — 37772 — 7772 + 3Z —7772 = 0. 

at at at at 

Thus the 2nd global invariant when passing (7772,7773) = (1,0) becomes 

yielding two separatrixes 



37773 - 777^ + 3Z 2 7772 = 2Z 3 , 



(5.7) 3t?7| = ra\- 3l 2 m 2 + 2Z 3 = (tt7 2 + 2Z)(777 2 - I) 2 . 

We note in passing that the two invariants ()5.4p .( i577|) are in fact the same spectral invariants 
we had before in (|2.3|) . which are now reformulated in terms of the traces 7772 and 7773. Thus, for 
example, the straightforward identity 

12m 4 - 77T7 2 . EE -4(4(J 2 - ft 2 - (?f - A8(b C ) 2 , 

reveals the relation betwen the trace-based invariant (|5.4p and the spectral invariant (|2.3b|) . 

In the phase plane (7772,7773), this consists of a closed curve for — 2Z < 7772 < Z and two open 
branches for 7772 > I. The phase plane analysis suggests that the global bounded solution exists 
if and only if the initial data satisfy (|5.5p and 

(7772, 7773) (0) G r, 



where 



Z - 7772 



r := I (m 2 , 777 3 )| |777 3 | < ^_ \fm 2 + 21, -21 < 1712 < I 



777 2 — Z / 1 

U |777 3 = — V m 2 + 21, 7772 > lj 



and the moving parameter Z is determined by (|5.5p . Also we can show that if initial data do not 
belong to T the solution becomes unbounded in finite time. 
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This figure depicts trajectories for system (15.60 in (7712, 777.3) plane, which contain the boundary 
of the non-trivial set T with (15. 5p . 

From above analysis it follows that the solutions remain bounded for all time if and only if 
the initial data Ao lies in the following set 

(5.8) A eu />0 {Snr}, 

where 

r :=riur 2 ur 3 , 

Ti := {(1712,7713) I \m 3 \ < =— v / 7Ti 2 + 21, -21 < m 2 < I}, 

v3 

T 2 := {(m 2 ,m 3 ) | |m 3 | = ^-y / m 2 + 21, -21 < m 2 < I}, 

m 2 - I 



T 3 := {(m 2 ,m 3 ) | m 3 = ^7712 + 21, m 2 > I}. 



We now turn to Theorem 11.31 Let r > be a moving parameter, we restrict attention to the 
parameterized surface m 2 + 21 = 3r 2 . Clearly the constraint m 2 > —21 is ensured for any real r. 
For any r > 0, the set S restricted on this surface is represented as 

4m 4 - 2m\ - 2?n 2 r 2 + 3r 4 = 0. 

This is a parabolic cylinder in the space (m 2 , 7713,7714). Applying the scaling property stated in 
Lemma 1.1, we may set r = 1, and denote the set S with constraint m 2 + 21 = 3 and mi = as 

4 

S := {A I 4m 4 - - 2m 2 + 3 = 0, ^1 = 0}, m k := ^ A^. 

fc=i 
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The first half of the set T|s is supported where —21 < m 2 < I, together with 7772 + 21 = 3, i.e., 
I = 3 ~™ 2 , leading to vri2 < 1. In this case, the restriction 

\m 3 \ <l - m 2 = m 2 = 1.5(1 - m 2 ), 

yields 

fil = {(m2,m 3 ) | 1 7713 1 < 1-5(1 - m 2 ), m 2 < l}. 
For reals A, m 2 > 0; the fact of no lower bound for m 2 suggests that any complex eigenvalue 
with zero divergence may well lie in S D Oi. 

The second r ^-constraint, supported on 777.3 = m 2 — I requiring 7772 > I = s ~™ 2 , i.e., 7772 > 1 
leading to 

2 = {(7772,7773) J 7773 = 1.5(7772 - 1), 777 2 > l}- 

The above set £ n [fii U £l 2 ] is 'fat'. Note the 3D case is similar to the special case / = which 
restricts to a large sub critical set. 



6. Appendix 

Finally we turn to interpretation of eigenvalues presented in Theorem 11.21 in terms of the 
subcritical sets in Theorem II. 3 1 or the equivalent set (|5,8D . 



{i} Two pairs of complex eigenvalues. The eigenvalues must be Ai = a + bi, X 2 = a — bi, A3 
—a + ci, A4 = —a — ci, where a, b, c S M. and b 7^ 0, c 7^ 0. A direct calculation gives 

m 2 = 4a 2 - 2b 2 - 2c 2 , 

7773 = 6a(c 2 — ft 2 ), 

7774 = 4a 4 + 2b 4 + 2c 4 - 12a 2 (6 2 + c 2 ), 

m \-l 2 = 3(7774 - ^777 2 ) = 3(6 2 - c 2 ) 2 - 24a 2 (& 2 + c 2 ), 

/2 = (4a 2 +6 2 + c 2 )2 + 126 2 c 2 > 

It follows that -21 < 777 2 < I and I > 4a 2 + b 2 + c 2 . Then 
I — m 2 



v / m 2 + 2/> 3(6 ~}L C ^ Vl2^ = 6\a\(b 2 + c 2 ) > |m 3 |. 



V3 V Z ' V3 
Thus we know S fl T = S f\ Ti. 

{ii} One pair of complex eigenvalues and two real eigenvalues. The four eigenvalues must be 
Ai = a + bi, X 2 = a — bi, A3 = —a + c, A4 = —a — c, where a, b, c G R and b ^ 0. Changing c in 
Case II to — ci we immediately obtain 

7T7 2 = 4a 2 - 26 2 + 2c 2 , 

777,3 = — 6a(6 2 + c 2 ), 

777 4 = 4a 4 + 26 4 + 2c 4 - 12a 2 (b 2 - c 2 ), 

ml -I 2 = 3(7774 - ^m 2 .) = 3(6 2 + c 2 ) 2 - 24a 2 {b 2 - c 2 ), 

I 2 = (4a 2 + b 2 - c 2 ) 2 - 126 2 c 2 . 

Suppose —21 < 7772 < ^ for / > 0, we distinguish two cases: 
(1) If 4a 2 + b 2 - c 2 > 0, then I < 4a 2 + b 2 - c 2 , and 

v^+2/ < 3(62 ~ c2) ^/l2a^ = 6)a(6 2 _ c 2 }| < | m;j |_ 
V 3 v 3 
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It becomes an equality if and only if c = 0. 

(2) If 4o 2 + b 2 - c 2 < 0, then I < c 2 - 4a 2 - b 2 , and I - m 2 = b 2 - c 2 - 8a 2 < 0. It's a contradiction 
with m2 < I- For 777,2 > I, the constraint 7713 = V m 2 + 21 leads to the relation (|6.ip . i.e., 

p = ^3m| - m| - 9m 2 (^n 4 - ^Tif.^ ^ + 108 - y7j™l^ = 0. 

Calculation shows that p = 432b 2 c 2 (b 2 + (2a - c) 2 ) 2 (b 2 + (2a + c) 2 ) 2 . So p = if and only if 
c = 0. Thus we know that S n T = S n T2, and two real eigenvalues must be equal. 



{hi} all the eigenvalues are real. Suppose the four real eigenvalues are a, b, c and —(a + b + c) 
a,b,c S R, then 7712 > 0. From the set S in (15.81) it follows 



^(mf, - I 2 ) = 7774 - ^7772 > 0, 

here we have used the inequality (a + (3 + 7 + 5) 2 < 4(a 2 + /? 2 + 7 2 + 5 2 ). These together lead 
to 777,2 > Thus if all the eigenvalues are real, then S n T = 5 n T 3 . 

Because of the homogeneousness, we can assume the four real eigenvalues are 1 + s, — l + w,— 1 
and 1 — s — w(if Aq 7^ 0). Let's do the following calculation. 

From 777,3 = !2 7^v / "i2 + 21 it follows that 

v3 

3ml = ( m 2 ~ lf(m 2 + 21) = m\- 3m 2 l 2 + 2l 3 , 
=> [3ml - m 2 + 3m 2 / 2 ] 2 = 4(Z 2 ) 3 . 



Using / 2 = — 3(777,4 — j^rn 2 ,) we have 



(6.1) p := ^3?77,1 — 7772 ~~ 9 m 2 — ^ + 1^8 ^777 4 — y7J m 2^ = 0- 

Calculation shows that 

p = _27(a + 2) 2 -u; 2 (s - w + 2) 2 (2s + w) 2 (s + 2w - 2) 2 (s + w - 2) 2 . 

So 77 = if and only if s = —2 orw = 0ors-w + 2= 0or2s + w = 0ors + 2w — 2 = 
or s + w — 2 = 0. They are all the same if we consider the homogeneousness and permutation. 
Now we know that the four eigenvalues must be in the form r(l + s, —1, —1, 1 — s). We claim 
that the range for s is [—2, 2]. 

(i) For -2 < s < 2, it's easy to check that Ae5nT 3 . 

(ii) For s > 2 or s < —2, we calculate Z, 777,2,7773 to obtain 



/ = s 2 -4, 777 2 = 4 + 2s 2 , 7773 = 6s 2 . 



So 



Hence 



7772 - I 



\Jm 2 + 21 



2 



7n 2 = ^(s 2 -4) 3 >0. 



V3 

We now can conclude that if all the eigenvalues are real, then S n T = S n 1^3 , and Ao must be 
in the form r(l + s, 1 — s — 1, — 1, ) (plus arbitrary permutation), where —2 < s < 2 and r > 0. 
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